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ABSTRACT 

The variation of the expansion rate of the Universe with time produces an evolution in 
the cosmological redshift of distant sources (for example quasar Lyman-a absorption 
lines), that might be directly observed by future ultra stable, high-resolution spectro- 
graphs (such as CODEX) coupled to extremely large telescopes (such as European 
Southern Observatory's Extremely Large Telescope, ELT). This would open a new 
window to explore the physical mechanism responsible for the current acceleration of 
the Universe. We investigate the evolution of cosmological redshift from a variety of 
dark energy models, and compare it with simulated data. We perform a Fisher matrix 
analysis and discuss the prospects for constraining the parameters of these models and 
for discriminating among competing candidates. We find that, because of parameter 
degeneracies, and of the inherent technical difficulties involved in this kind of observa- 
tions, the uncertainties on parameter reconstruction can be rather large unless strong 
external priors are assumed. However, the method could be a valuable complementary 
cosmological tool, and give important insights on the dynamics of dark energy, not 
obtainable using other probes. 

Key words: Cosmology: theory - cosmological parameters - Cosmology: observations 
- Galaxies: distances and redshifts 



1 INTRODUCTION 

' The discovery of the current accelerated expansion of the 
, Universe (Riess et al. 1998; Perlmutter et al. 1999) has 
prompted a renewed interest towards classic cosmological 
tests. The measurement of the expansion rate of the Uni- 
verse at different redshifts is crucial to investigate the cause 
of the accelerated expansion, and to discriminate candidate 
models. Until now, a number of cosmological tools have been 
successfully used to probe the expansion and the geome- 
try of the Universe. The position of acoustic peaks in the 
cosmic microwave background (CMB) angular power spec- 
trum provides a powerful geometrical test, showing that the 
space curvature of the Priedmann-Robertson- Walker (FRW) 
metric is nearly flat (Spergel et al. 2006). A similar test is 
performed through the detection of baryon acoustic oscil- 
lations (BAG) in the power spectrum of matter extracted 
from galaxy catalogues. The luminosity distance of type la 
supernovae and other standard candles allows to constrain 
the value of the expansion rate at different redshifts, typi- 
cally up to z ~ 2 (Riess et al. 2006). 
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Currently, however, very little is known about the de- 
tailed dynamics of the expansion. Depending on the underly- 
ing cosmological model, one expects the redshift of any given 
object to exhibit a specific variation in time. An interesting 
issue, then, is to study whether the observation of this vari- 
ation, performed over a given time interval, could provide 
useful information on the physical mechanism responsible 
for the acceleration, and be able to constrain speciflc mod- 
els. This is the main goal of this paper. In addition to being 
a direct probe of the dynamics of the expansion, the method 
has the advantage of not relying on a determination of the 
absolute luminosity of the observed source, but only on the 
identification of stable spectral lines, thus reducing the un- 
certainties from systematic or evolutionary effects. The pos- 
sible application of this kind of observation as a cosmological 
tool was flrst pointed out by Sandage (1962). However, the 
tininess of the expected variation (typically, a shift of less 
than a cm/s over a period of a year) was deemed impos- 
sible to observe at the time. The importance of this test 
was stressed again over the past decades by other authors 
(e.g. Lake 1981, Riidiger 1980); more recently, Loeb (1998) 
re-assessed its feasibility and concluded that, given the ad- 
vancement in technology occurred over the last forty years, 
it is conceivable to expect that a measurement of the redshift 
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variation in the spectra of some suitable source (most no- 
tably the Lyman-a absorption lines of distant QSOs) could 
be detected in the not too distant future. Recently, Lake 
(2007) showed that measuring the time evolution of redshift 
would be a way to check the internal consistency of the un- 
derlying cosmological model, and to map the equation of 
state of dark energy. 

With the foreseen development of extremely large obser- 
vatories, such as the European ELT, with diameters in the 
range 30-100 in, and the availability of ultra-stable, high- 
resolution spectrographs, the perspectives for the observa- 
tion of redshift variations look very promising. For example, 
Pasquini ct al. (2005, 2006) pointed out that the CODEX 
(COsmic Dynamics Experiment) spectrograph should have 
the right accuracy to detect the expected signal by monitor- 
ing the shift of Lyman-a lines of distant {z > 2) QSOs over 
a period of some decades. 

A previous investigation of the expected cosmological 
constraints from this kind of observations was performed 
by Corasaniti et al. (2007). That work, however, only made 
predictions for a very restrictive set of models and assump- 
tions: firstly, it only explored the case when the dark en- 
ergy is a standard cosmological constant (i.e. a component 
with constant equation of state w = p/p — —1) plus a few 
non-standard models (the Chaplygin gas, and a model that, 
although named interacting daxk energy, in fact only gen- 
eralises a constant dark matter equation of state, leaving 
the cosmological constant untouched); secondly, as it will 
be shown in more detail later, the observational strategy 
considered in that paper seems rather optimistic. 

In this paper, we aim to give an overview of the theoret- 
ical predictions for the most popular and still viable models 
introduced in the recent past to explain the observed accel- 
erated expansion. These models either fall within the broad 
category of "dark energy" , introducing an unknown, smooth 
and gravitationally repulsive component, or invoke a mod- 
ification of the theory of gravity (for a review of possible 
explanation for the accelerated expansion see, e.g., Peebles 
& Ratra (2003)). Our first goal is to investigate whether the 
dynamics of any of these standard and non-standard dark 
energy models shows interesting features that could be con- 
strained by future observations of the redshift variation. We 
do not restrict ourselves a priori to models with fiat geom- 
etry, and we fully take into account correlations between 
parameters when assessing the expected uncertainties. 

This paper is organised as follows. In Sect. 2 we review 
the basic equations that describe the redshift variation with 
time of a source, in an expanding universe. We comment on 
the possibility of detecting the effect in the near future in 
Sect. 3. In Sect. 4 we study the expected apparent veloc- 
ity shift for a rmmber of non-standard dark energy models. 
Finally, we discuss the main results and conclusions of our 
work in Sect. 5 and 6. 



2 BASIC THEORY 

We start by reviewing the basic theory necessary to de- 
rive the expected redshift variation in a given cosmologi- 
cal model. We assume that the metric of the Universe is 
described by the FRW metric. The observed redshift of a 



given source, which emitted its light at a time ts, is, today 
(i.e. at time to), 
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The observed redshift variation of the source is, then, 
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which can be re-expressed, after an expansion at first order 
in At/t, as: 



Azs ~ Ato 
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Clearly, the observable Azs is directly related to a change 
in the expansion rate during the evolution of the Uni- 
verse, i.e. to its acceleration or deceleration, and it is then 
a direct probe of the dynamics of the expansion. It van- 
ishes if the Universe is coasting during a given time in- 
terval (i.e. neither accelerating nor decelerating). We can 
rewrite the last expression in terms of the Hubble parame- 
ter H{z) = d{z)/a{z): 



Azs = HoAto l + zs- 



H(zs 
Ho 



(5) 



The function H{z) contains all the details of the cosmologi- 
cal model under investigation. Finally, the redshift variation 

can also be expressed in terms of an apparent velocity shift 
of the source, Av — cAzs/{l + Zs). 



3 CAN THE VELOCITY SHIFT BE 
OBSERVED? 

The latest studies on the feasibility of detecting a time evolu- 
tion of the redshift are those by Pasquini et al. (2005, 2006). 
The most promising approach is to look at the spectra of 
Lyman-Q forest absorption lines, which are very stable and 
basically immune from peculiar motions. This is a scenario 
that might be achieved in the next decades, when extremely 
large telescope (such as the ELT) will collect a large num- 
ber of photons, and high-resolution spectrographs (such as 
CODEX) will be able to measure tiny shifts in spectral lines 
over a reasonable amount of time, typically of order of few 
decades. 

According to Monte Carlo simulations discussed by 
Pasquini et al. (2005, 2006), the accuracy of the spectro- 
scopic velocity shift measurements expected by CODEX can 
be modelled as: 
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cm/s; (6) 



here, S/N is the signal to noise ratio for pixels of 0.0125 A, 
Nqso is the inimbcr of QSOs spectra observed and zqso is 
their redshift. Based on the currently known QSOs brighter 
than m ^ 16.5 in the redshift range 2 < « < 4, Pasquini 
et al. (2005, 2006) assumed to observe either 40 QSOs with 
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S/N ratio of 2000 each, or 30 QSOs with S/N of 3000, re- 
spectively. In fact, again according to Pasquini et al. (2005), 
a CODEX-hke experiment, coupled to a 60 m telescope with 
approximately 20% total efficiency, would give a cumulative 
S/N of 12000 for a single QSO, requiring roughly 125 hours 
of observation to get a S/N of 3000 on that spectrum. Then, 
starting with 10 hours of observation per night, and taking 
into account a 20% use of the telescope, and a 90% of actu- 
ally usable data, one finds that 40 spectra can be measured 
with that S/N in roughly 7.6 years (this time would actu- 
ally increase to about 15 years if the telescope aperture is 
smaller, e.g., 42 m instead of 60 m). Then, it seems that a 
reasonable time span to perform a second observation of the 
same spectra might be 30 years. 

Based on the previous considerations, our study will 
then be conducted assuming a future dataset containing 
a total of 40 QSOs spectra (uniformly distributed over 5 
equally spaced redshift bins in the redshift range 2-5), with 
a S/N = 3000, observed twice over a time interval of 30 
years. We note that a previous study performed by Corasan- 
iti et al. (2007) seems to make optimistic claims, i.e. that 240 
pairs QSO's spectra can be observed with S/N = 3000 over 
a time span of 10 years. It is easily shown (with the same 
arguments as above) that with these figures the required ob- 
serving time would actually be roughly 90 years for each of 
the two epochs. Moreover, even assuming an increase of the 
number of QSOs in the redshift range 2-4 with future large 
all-sky photometric surveys, it seems quite difficult to pre- 
dict an order of magnitude increase over the current known 
objects, which axe 25 (Pasquini et al. 2005). 

Using the expected error bars from Eq. 6, we can pre- 
dict the level of accuracy that can be reached in the recon- 
struction of the parameters for a comprehensive set of dark 
energy models. Furthermore, we can predict whether any of 
these models can be distinguished from the standard ACDM 
scenario. We will now look into these problems. 



less stated otherwise, throughout our paper we assume for 
each class of models the best fit values found in that work, 
and vary the parameters within their 2a uncertainties. 

All the predictions on the time evolution of redshift pre- 
sented in the following were derived assuming observations 
performed over a time interval Ato = 30 years. From Eq. 5 
it is clear that the expected velocity shift signal increases 
linearly with Ato, so that it is straightforward to calculate 
the expected signal when a different period of observation is 
assumed. Fig. 1 shows our predictions for the cosmological 
models discussed in the following, along with simulated data 
points and error bars derived from Eq. 6. 



4.1 The standard cosmological model (ACDM) 

We start our analysis by first setting out the predictions for 
the current standard cosmological model. In the simplest 
scenario, the dark energy is simply a cosmological constant, 
A, i.e. a component with constant equation of state w — 
p/p = —1. Flatness of the FRW metric is usually assumed, 
but in general one can parametrize the deviation from the 
zero-curvature case in terms of the parameter fife s 1 — ft 
where Q is the total density of the Universe in units of the 
critical value pc = 3Hq /8itG. 

The Hubble parameter evolves according to the Pried- 
mann equation, which, for this model, is: 



(7) 



where Q.m and SIa parameterize the density of matter and 
cosmological constant, respectively. When flatness is as- 
sumed, n — Qm + = 1, and the model has only one free 
parameter, Q,a- The current best fit value from cosmological 
observations is Oa ~ 0.73 ±0.04 in the flat case (Davis et al. 
2007), while relaxing the assumption of flatness results in a 
preference for slightly closed models, with Qfc = — O.OliO.Ol 
(Spergel et al. 2006). 



4 PREDICTIONS FOR DARK ENERGY 

MODELS 

All of the models we will consider in this paper are currently 
viable candidates to explain the observed acceleration, i.e. 
they have not been falsified by available tests of the back- 
ground cosmology. Clearly, some models may be preferred 
with respect to others based on some statistical assessment 
of their "economy" , i.e. the fact that they fit the data well 
with a smaller number of parameters. Given the current sta- 
tus of cosmological observations, there is no strong reason 
to go beyond the simple, standard cosmological model with 
zero curvature and a cosmological constant A (except for the 
conceptual problems arising when one attempts to reconcile 
its observed value with some estimate derived from funda- 
mental arguments, sec, e.g., Weinberg 1989). For the scope 
of the present paper, however, it is interesting to explore as 
many models as possible, since future observations of the 
time variation of redshift could reach a level of accuracy 
which could allow to better discriminate competing candi- 
dates, and to understand the physical mechanism driving 
the expansion. We refer the interested reader to the paper 
by Davis et al. (2007), which discusses the constraints on 
most of the models that we will focus on in our study. Un- 



4.2 Dark energy with constant equation of state 

The next step is to allow for deviations from the simple 
w = —1 case, introducing a component with an arbitrary, 
constant value for the equation of state. The accelerated 
expansion is obtained when w < —1/3. The Hubble param- 
eter for this generic dark energy component with density 
fide then becomes: 



de 



-3{l-|-n;) ■ 



(8) 



The currently preferred values of w in this models still in- 
clude the cosmological constant case: w = —1.01 ± 0.15 
(Davis et al. 2007). 



4.3 Dark energy with variable equation of state 

If the equation of state of dark energy is allowed to vary 
with time, then the Hubble parameter is: 



n. 



de 



i J da(l-\-w(a)) / a 



(9) 
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In this case, one has to choose a suitable functional form for 
w(o), which in general involves a parametrization. The most 
commonly used (Chevallier & Polarski 2001; Linder 2003) is: 

■w{a) = Wo + Wail — a), (10) 

although different approaches can be used. Clearly, the ex- 
act form of w{a) with time will lead to completely different 
evolution for the dark energy component. 

4.4 Interacting dark energy 

In interacting dark energy models the dark components in- 
teract through an energy exchange term. The conservation 
equations for matter and dark energy can be written in a 
very general way as 

p'm + SHpm = SHpm, (11) 
Pde + 3Hpde{l+w)=-5Hpm, (12) 

so that the total energy-momentum tensor is conserved. 
Whenever 5 is a non-zero function of the scale factor, the 
interaction causes pm and pde to deviate from the standard 
scaling, and the mass of the particles is not conserved. This 
non-standard behaviour has been parametrized (Dalai et al. 
2001; Majcrotto ct al. 2004) by the relation pde/pm = Aa^, 
where A = f2de/(l — ^de — ^k) and the density parameters 
are the present quantities. The Hubble parameter or this 
model then reads 

which reduces to the uncoupled case for ^ = — 3w. This 
model also includes all late-time scaling solutions. We also 
note that this model is a genuinely interacting dark energy, 
unlike the one discussed in Corasaniti ct al. (2007) which is 
actually a generalised dark matter (thus more similar to the 
model we discuss in 4.8). 

4.5 DGP models 

The Dvali-Gabadadze-Porrati (DGP) model (Dvali et al. 
2000) provides a mechanism for accelerated expansion which 

is alternative to the common repulsive-gravity fluid ap- 
proach: within the context of brane-world scenarios, the 
leaking of gravity in the bulk, above a certain cosmologi- 
cally relevant physical scale, is responsible for the increase 
in the expansion rate with time. The only parameter of this 
class of models is r, the length at which the leaking occurs, 
which defines an adimensional parameter Qr = l/(4r^i/o)- 
The Hubble parameter then reads: 




where = 1 — fi/t — 2\/f2rVl ~ ^k- 



4.6 Cardassian models 

Another possibility originating from the brane-world sce- 
nario is that of a so-called Cardassian expansion (Freese & 
Lewis 2002) resulting from a modification of the Priedmann 
equation, with the introduction of a term that depends non 



linearly on the average density of the Universe, assumed 
to be composed only of matter. This additional term, phe- 
nomenologically, is equivalent to the introduction of a dark 
energy component, with a scaling law oc a~^", where n is 
completely equivalent to the quantity w -|- 1 of the dark en- 
ergy models with constant equation of state. More inter- 
esting are the so-called "modified polytropic Cardassian" 
models, which have an extra parameter q and a Ftiedmann 
equation: 




4.7 Chaplygin gas 

There are a few models which attempts to explain both 
structure formation and the current acceleration of the Uni- 
verse with a single "dark fluid", whereas the standard sce- 
nario relies on two separate dark contributions to the stress- 
energy tensor (a dark matter and a dark energy component). 
A well-studied case is the so called Chaplygin gas (Kamen- 
shchik et al. 2001), whore the unified dark component has 
equation of state p = Ap~'^ with A > 0. 

The expansion rate in this model is governed by the 
equation: 

(S)'^^^a-)(-S^)"""' a., 

with the definition A — A/p^f'^"', where po is the total 
density of the Universe at the present. The so-called "stan- 
dard case" for the Chaplygin gas is obtained for the choice 
7 = 1 (which, however, is not a good fit to current data), 
while for 7 = the model recovers the standard cosmological 
constant case with fim, = {1 — iik){i — A). 



4.8 AfHne equation of state 

An interesting possibility to consider is that the dark en- 
ergy is modelled by a generic, barotropic equation of state 
P = P{p)i as discussed in Chiba et al. (1997); Visser (2004); 
Ananda & Bruni (2006). In particular, the case where the 
Taylor expansion of an arbitrary equation of state of that 
sort is truncated to first order, e.g. p = po + ctp, has recently 
been investigated by Balbi et al. (2007), who also derived 
cosmological constraints on the parameters of the model. It 
is interesting to note that such an afiine equation of state 
can be used to describe a simple unified dark matter model, 
with a time evolution of the background density given by 
p{a) = PA + {po - pa)o~*^^+"^ , where pA = -po/(l + a) and 
Po is the dark energy fluid energy density at present. The 
Hubble parameter is given by: 

where = {po — Pa)/ Pc', for a = this model recovers the 
standard ACDM case. 
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5 RESULTS AND DISCUSSION 

We used the expected error bars from Eq. 6 to perform a 
Fisher matrix analysis (Tegmark et al. 1997), leading to an 
estimate of the best possible constraints on the parameters 
of dark energy models. All our predictions are based on 30 
years of observation, assuming that the fiducial values for 
the parameters of each dark energy model are those which 
best fit current observations (as from Davis et al. 2007, Balbi 
et al. 2007, Majerotto et al. 2004). 

The Fisher matrix formalism allows one to estimate the 
best possible accuracy attainable on the determination of 

the parameters of a certain model. Specifically, given a set 
of cosmological parameters pi, i = l,...,n, and the corre- 
sponding Fisher matrix Fij (that is easily calculated based 

on a theoretical fiducial model and the assumed data errors 
from Eq. 6), the best possible la error on pi is given by 

1/2 

Api = Cj/ , where the covariancc matrix dj is simply the 
inverse Fisher matrix: Ca = It is a well-known fact 

(see, e.g. Bond et al. 1997) that, when estimating the ex- 
pected errors Api on each parameters, one has to be careful 
about existing correlations with the other parameters. Dur- 
ing the inversion process, non-null off-diagonal Fisher matrix 
elements propagate in the diagonal elements of the covari- 
ance matrix, giving their contribution to the estimated un- 
certainties. Neglecting this contribution, for example by sim- 
ply taking the inverse of the elements of the Fisher matrix (a 
process which is equivalent to assume perfect knowledge of 
all the other parameters) usually results in severely under- 
estimated errors. The existence of correlations among esti- 
mated parameters is a manifestation of degeneracies: when 
more and more parameters are allowed to vary, and they 
have similar effects on the observable quantities, it becomes 
increasingly difficult to constrain each parameter. 

From the above considerations, two crucial conse- 
quences arise: first, dark energy models with intrinsically 
less free parameters will have an advantage with respect to 
models with more free parameters; second, not taking prop- 
erly into account the degeneracies among parameters (for 
example, by assuming perfect knowledge on some of them) 
will lead to wrong estimates of the errors. For this reason, we 
allowed all the parameters which are specific of a given dark 
energy model to vary simultaneously in our analysis (i.e. 
we inverted the full Fisher matrix when estimating errors). 
However, since it would be pointless not to assume any prior 
knowledge, we fixed the parameters not directly related to 
dark energy (such as, for example, the baryon density) to 
their current best fit value. 

As an example of what we just discussed, let us consider 
first the standard ACDM case. Assuming that the fiducial 
model has = 0.7 and 0^ = and that both and fi/c 
can vary, we find AQa = 0.2 and Afifc = 0.25 at Icr. When 
we fix Qk ~ 0, the bound on Qa becomes 0.007 (at la). 
If dark energy is modelled by a constant equation of state 
(with a fiducial value w = — 1) and the flatness constraint is 
imposed we find a looser bound on the dark energy density, 
A^lde = 0.016, and quite a large error on the equation of 
state, Aw = 0.58. This clearly shows that different assump- 
tions on the knowledge of any parameter has an infiuence 
on all the others. Having noted this, however, we also note 
that the parameter flfc can be much better constrained using 



external datasets, such as the CMB anisotropy. For the sake 
of simplicity, then, we will assume flatness in the following. 

The prospect of detecting departures from the standard 
ACDM case could in principle be one of the real assets of 
observing the time evolution of redshift, and is thus worth 
closer investigation. Since the simulated data used in our 
analysis assume that QSOs are used as a tracer of the red- 
shift evolution, we expect that the more constrained models 
will be those that have the largest variability in the redshift 
range 2 < 2: < 5. Of course, additional information could be 
obtained using other observational windows. We will com- 
ment on this issue later on. 

The DGP model is the one for which we obtained the 
tightest constraints in our analysis: Aflr = 0.0027 at Icr, 
assuming Qr = 0.13 as a flducial value. This is not only 
due to the strong dependence of the velocity shift on fir 
(see Fig. 1), but also to the simplicity of the model, which 
depends on only one parameter (in this respect, this is the 
simplest model, together with the standard flat ACDM). In 
general, as we already discussed, it is to be expected that 
models with less parameters perform better. 

The issue of degeneracies is crucial when looking into 
the behaviour of models having more than one free param- 
eters. We start by discussing the Cliaplygin model. When 
both A and 7 vary freely, no interesting constraint can be 
obtained observing the velocity shift with the assumed QSO 
data: we find = 0.42 and A7 = 1.4 (for the fiducial val- 
ues A = 0.7 and 7 — 0.2). Fixing A, on the contrary, results 
in a very tight bound on 7: A7 = 0.008. 

We then address the interacting dark energy and the 
affine equation of state. As it is clear from Fig. 1 these mod- 
els show a large variability in the redshift range we are ex- 
ploring: this is to be expected, since in both models the 
matter-like component departs from the usual scaling, 
giving a distinctive signature when one looks at higher and 
higher redshifts. We could then naively expect tight bounds 
on the parameters of these models from the observation of 
velocity shift. Unfortunately, again, this is not the case, due 
to strong degeneracies among the parameters. For example, 
when S^A is allowed to vary, for the affine equation of state 
we find Aa = 0.05 (for a fiducial value a = 0.02). How- 
ever, if we suppose that one of the parameters is known 
from other observations, the constraints on the others nar- 
row considerably. If Qa is known, the affine parameter a 
can be reconstructed with an error Aa — 0.005. If, in addi- 
tion, we also know that w = —1 we find A^ = 0.06 for the 
interacting dark energy model (for a fiducial value ^ = 3). 

The other models do not seem to have very interest- 
ing signatures to be exploited, at least in the redshift range 
considered in our analysis. The worst bounds are expected 
for the Cardassian models, which has not a large parameter 
dependence in the redshift range 2 < « < 5 (see Fig. 1). 
In the case of a varying equation of state, the simple pa- 
rameterization adopted in our study shows significant de- 
viations from the fiducial model only when rather extreme 
values of Wa (already excluded by current observations) are 
assumed. However, alternative parametrizations might lead 
to different scalings, and dark energy might dominate in the 
redshift window probed by Lyman-Q forests. Models with 
this behaviour can in principle exist (see, for example, Do- 
delson et al. 2000, Frieman et al. 1995, Caldwell & Linder 
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2005, Linder 2006). A recent review on the dynamics of dark 
energy models is Linder (2007). 

A different aspect that we dealt with in our study was 
the potential of this cosmological tool to increase the ability 
to discriminate different models, when combined to other ob- 
servations. Fig. 2 shows the comparison among the predicted 
velocity shifts for all the models described earlier, assuming 
parameter values that are a good fit to current cosmolog- 
ical observations (including the peak position of the CMB 
anisotropy spectrum, the SNe la luminosity distance, and 
the baryon acoustic oscillations in the matter power spec- 
trum). In other words, the models shown in Fig. 2 cannot 
be easily discriminated using current cosmological tests of 
the background expansion. If we assume that the ACDM 
model is the correct one, and simulate the corresponding 
data points for the velocity shift, using a test we can 
quantify how well we can exclude the competing models 
based on their expected signal. As it is clear from Fig. 2, 
some models can be excluded with a high confidence level. 
In particular, we find that the Chaplygin gas model and the 
interacting dark energy model would be excluded at more 
than 99% confidence level, and that the afHne model would 
be out of the 1 a region. 

As a final caveat, we note that our results were obtained 
assuming an equal number (8) of QSOs for each of 5 red- 
shift bins in the range z =2-5. (Such a uniform distribution 
was also assumed in the simulations performed by Pasquini 
et al. 2006.) This explains the fact that our error bars sig- 
nificantly decrease with redshift (see Fig. 1). Assuming a 
decreasing number of QSOs at higher redshifts (undoubt- 
edly, a somewhat more realistic assumption) would result in 
a slight increase in the error bars for those bins. For exam- 
pic, assuming 3 QSOs instead of 8 in the highest redshift 
bin would increase the error bar for that bin of a factor 1.6. 
However, since the largest variability in theoretical predic- 
tions is precisely at high redhisfts, we would not expect our 
main conclusions to change dramatically. 



difference too small to produce any interesting constraint on 
cosmological parameters. 

Despite its inherent difficulties, the method has many 
interesting advantages. One is that it is a direct probe of 
the dynamics of the expansion, while other tools (e.g. those 
based on the luminosity distance) are essentially geometrical 
in nature. This could shed some light on the physical mecha^ 
nism driving the acceleration. For example, even if the accu- 
racy of future measurements will turn out to be insufficient 
to discriminate among specific models, this test would be 
still valuable as a tool to support the accelerated expansion 
in an independent way, or to check the dynamical behaviour 
of the expansion expected in general relativity compared to 
alternative scenarios. Furtliermore, despite being obscrva- 
tionally challenging, the method is conceptually extremely 
simple. For example, it does not rely on the calibration of 
standard candles (as it is the case of type la SNe) or on a 
standard ruler which originates from the growth of perturba- 
tions (such as the acoustic scale for the CMB) or on effects 
that depend on the clustering of matter (except on scales 
where peculiar accelerations start to play a significant role). 
Moreover, the errors on the measured data points decrease 
linearly with time and can become significantly smaller over 
only a few decades of observations. Finally, it is at least con- 
ceivable that suitable sources at lower redshifts than those 
considered in this work could be used to monitor the veloc- 
ity shift in the future. This would be extremely valuable, 
since some non-standard models have a stronger parameter 
dependence at low and intermediate redshifts (see Fig. 1), 
that could be exploited as a discriminating tool. In Loeb 
(1998) speculative possibilities of using other sources have 
been indicated, like masers in galactic nuclei, extragalactic 
pulsars or gravitatioiially lensed galaxy surveys: this would 
further extend the lever arm in redshift space and increase 
the ability of constraining models. Exploring the feasibility 
of such proposals is beyond the scope of this paper, but it 
may certainly be an interesting topic for further studies from 
observers. 



6 CONCLUSIONS 

In summary, we found that the measurement of the veloc- 
ity shift with future extremely large telescopes and high- 
resolution spectrographs could provide interesting informa- 
tion on the source of cosmic acceleration, which would com- 
plement other, more traditional cosmological tools. From our 
analysis, it is also clear that the observation of velocity shift 
alone can be affected by strong parameter degeneracies, lim- 
iting its ability to constrain cosmological models. Contrarily 
to other analyses, then, our conclusion is that the uncer- 
tainties on parameter reconstruction (particularly for non- 
standard dark energy models with many parameters) can 
be rather large unless strong external priors are assumed. 
When combined with external inputs, however, the time evo- 
lution of redshift could discriminate among otherwise indis- 
tinguishable models. 

We also found that a reasonable time span to perform 
the comparison on the redshift evolution seems to be roughly 
three decades. Shorter time spans seem unrealistic, given 
the time needed to observe QSOs spectra with the neces- 
sary S/N. A shorter interval would also make the redshift 
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Figure 1. The apparent spectroscopic velocity shift over a period Ato = 30 years, for a source at redshift Zs, for the models described 
in the text. Prom top to bottom and from left to right: the ACDM (Sect. 4.1), dark energy with a constant equation of state (Sect. 
4.2), dark energy with varying equation of state (Sect. 4.3), interacting dark energy (Sect. 4.4), DGP (Sect. 4.5), Cardassian (Sect. 4.6), 

generalized Chaplygin gas (Sect. 4.7), dark fluid with an afflne equation of state (Sect. 4.8). All the other parameters are fixed at their 
best fit value. The decreasing error bars are due to the assumption of a uniform distribution of QSOs over the entire redshift range (see 
Sect. 5 for a discussion). 
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Figure 2. The predicted velocity shift for the models explored in this work, compared to simulated data as expected from the CODEX 
experiment. The simulated data points and error bars are estimated from Eq. 6, assuming as a fiducial model the standard ACDM model. 
The other curves are obtained assuming, for each non-standard dark energy model, the parameters which best fit current cosmological 
observations. 



